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Abstract 



We report on a Monte-Carlo study of two-dimensional Ginzburg-Landau 
superconductors in a magnetic field which finds clear evidence for a first- 
order phase transition characterized by broken translational symmetry of the 
superfluid density. A key aspect of our study is the introduction of a quantity 
proportional to the Fourier transform of the superfluid density which can 
be sampled efficiently in Landau gauge Monte-Carlo simulations and which 
satisfies a useful sum rule. We estimate the latent heat per vortex of the 

melting transition to be ~ 0.38A:bTa// where Tm is the melting temperature. 
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In mean-field theory type II Ginzburg-Landau superconductors in a magnetic field have 
an unusual second-order phase transition. In the low-temperature (T < T^^^) phase dis- 
covered by Abrikosov the zeros of the superconducting order parameter (vortices) form 
a lattice and the system exhibits both broken translational symmetry and off-diagonal long- 
range order (ODLRO). Unusual aspects of the transition are related to the Cooper-pair 
Landau level structure which causes the mean field instability of the disordered phase to 
occur simultaneously at T^^^ in a macroscopic number of channels. However, the nature of 
this phase transition is qualitatively altered by thermal fluctuations. Interest in the effect 
of thermal fluctuations on the thermodynamic properties of type II superconductors has 
increased since the discovery of high-temperature superconductors which have an unusu- 
ally short coherence length so that fluctuation effects are important over a relatively wide 
temperature interval surrounding T^^^ . 

For D dimensional superconductors fluctuations in a magnetic field at temperatures well 
above T^^^, where different channels are independent, are like those of a -D — 2 dimen- 
sional system at zero magnetic field suggesting that the mean-field phase transition to 
the Abrikosov state will be destroyed by fluctuations for D < 4. High-temperature pertur- 
bative expansions , even when evaluated to high-order where coupling between different 
channels becomes important, show no evidence of a transition for D = 3 or D = 2 between 
the high-temperature fluid state and Abrikosov's vortex-lattice state. The results of Monte- 
Carlo simulations for D = 2 have been controversial. Tesanovic and Xing and Kato and 
Nagaosa find evidence for a phase transition at a temperature below T^^^ while O'Neill 
and Moore have concluded that the Abrikosov phase transition is suppressed by thermal 
fluctuations. In this letter we present the results of a Monte-Carlo simulation for D = 2 
in which we find unambiguous evidence for a first-order phase transition. 

The free energy density of a Ginzburg-Landau superconductor is given by 

fm = airm' + + ^l(-^^V - 2eA)^\'. (1) 

[F = J (i^r/[^(r)].) The quadratic terms in Eq. ([l|) are minimized by order-parameters 
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which correspond to a lowest Landau level (LLL) wavefunction for the Cooper pairs. It 
follows that the mean- field theory superconducting instabihty occurs at T^^^ (aniT^^^) = 
0; an = a + heB /m* ) for all Cooper pair states which are in the LLL but only at much lower 
temperatures for channels corresponding to higher Landau level Cooper pair wavef unctions. 
In this work we adopt the LLL approximation in which we assume that fluctuations in 
higher Landau level channels can be neglected [^|lOl and consider only the two dimensional 
limit where variations of the order parameter along the z direction can be neglected. In 
the LLL approximation the order parameter is defined up to an overall scale factor by its 
zeros, i.e. by the positions of the vortices. (This property has been used by Tesanovic and 
collaborators to develop many useful insights.) This limit applies to films thinner than 
a coherence length and to layered systems when the inter-layer coupling can be neglected. We 
choose the Landau gauge {A = (0, Bx, 0)) and apply quasi-periodic boundary conditions to 
the order parameter inside a finite system with lengths and Ly. (For thin films, especially 
those formed of strongly type II materials it is a good approximation to ignore fluctuations 
in the vector potential A) The order parameter \l/(r^ can then be expanded in the form, 

m = ( '"^'"^'^^ )^^^E^.[E(-^.^.)'^^^K)'^^^exp(z^X,,/f^) exp(-(x - X,,r/Af)] 

(2) 

In Eq. (^ Xj^s = j2TTi'^/Ly + sLx, (? = hc/2eB, s runs over all integers and j runs from 1 
to A^^ = LxLy/ {2ttP) which must be chosen to be an integer. 

A central role in our study is played by the superfluid-density spatial correlation function, 
whose Fourier transform is defined by 

XsFDik) = / d'r ^d'f'mr)\'\iJ{f'\')exp[^k■ (f-r^)] (3) 

LxLy J J 

We evaluate XsFoik) by expressing it in terms of 

^(^) = ^ E Q.Q.^..-.i-n, eM-^UXn + X,,)/2\ (4) 

where for a finite system k = 2'n'{nx/Lx,ny/Ly) , 6j = 1 if j is a multiple of A^^ and 
is zero otherwise, and Xj = Xj^q. (Note that Aq = A{k = 0) is proportional to the 
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integrated superfluid density.) A(fc) is conveniently sampled in our Landau gauge Monte 
Carlo simulations and 

Xsrnik) = ^(^^^)^exp[-fc^£V2](|A(fc)r) (5) 

— * 

Moreover A(fc) satisfies the following sum rule for each configuration of the Ginzburg-Landau 

system, 

-^Y.mk)\'-im = o (6) 

where A(A;) = A(A;)/Ao. Note that A(A;) depends only on the distribution of |\Ef(f^p and 
not on its overall magnitude. Eq. reflects the LLL restrictions on the superfluid density 
distribution. (For a finite system, both Ux and in the sum over k in Eq. (|^) range over 
any consecutive values.) 

In the vortex-liquid state XsFoik) should be a smooth function of wavevector and if the 
sum over k in Eq. is to converge we must have lim|fc|^oo |A(A;)p A^^^ and hence that 

, hm xsFDik) = exp[-A;2£V2] (7) 



It is readily verified |jTT|] that Eq. (0) is satisfied for all /c 7^ when T ^ T*^^ and the vortex 

fluid is completely uncorrelated. On the other hand, in a vortex-lattice state A(A;) = Ao^g^j 

where G is a reciprocal lattice vector. To see that Eq. is satisfied in this case note that 

there are A^,^ wavevectors per BriouUin zone in the Abrikosov state. Eq. (^) tells us that 

averages of |A(/i;)p over large areas of reciprocal space yield A^^^ irrespective of the degree 

of correlation among the vortices. A^(^(|A(/c)p) — 1 provides a very convenient measure of 

the degree of vortex correlation in a system. 

We can express the Ginzburg-Landau free energy in terms of |A(A;)p as follows 

m\ Epil^ , /3[A]Ag ^ 

"'^ = — — = [sgn(aH)Ao H J (8) 



ksT ksT 



where 



l3[A]^j:\m\'exp[-^]. (9) 



and g = ani'^^'^Lz/ I3KbTY/'^ . /9[A] has its minimum value in the Abrikosov state and 
increases as the vortex positions become less correlated. It is readily verified that in 
the uncorrelated vortex fiuid /3[A] = 2 while for the triangular lattice Abrikosov state 
/9[A] = ~ 1.159595. (This relatively weak variation in (3 was exploited recently [0 by 
Tesanovic et ai.) We regard /3[A] and Aq as the two intensive thermodynamic variables 
which characterize the state of the LLL Ginzburg-Landau system. We can define an en- 
tropy which measures the function-space volume associated with a given /3[A] and Aq by 
Sf,{l3,Ao) = kBHW{(3, Ao)) where 

W{P, Ao) ^ ^\^M}^llhf^ Y[ J dC,dC,5{(3 - m)5{A, - Cfi,) (10) 

With this definition the free energy, (3 and Aq at any value of g can be determined by 
minimizing 

Fp{P, Ao) = Ep{/3, Ao) - TS{/3, Ao) (11) 

with respect to (3 and Ao. {Fp is extensive so fiuctuations become negligible in the ther- 
modynamic limit.) We will use Eq. (|TI|) to interpret the Monte Carlo results discussed 
below. 

Using the Metropolis algorithm we have determined distribution functions for several 
quantities including Ep, Aq, i3[A], and |A(A;)p as a function of both g and N^. Finite 



system shapes have been chosen to accommodate perfect triangular lattices. For all simu- 
lations the order parameter was initialized to the Abrikosov lattice value and the first 10'^ 
Monte Carlo steps were discarded. Some typical results for (|A(A;)p) at T < T*^^ are shown 
in Fig. (|l]). At g^ = 30 the vortex fiuid has developed strong correlations. For A^,^ = 120, 
N^{\A{G)\'^) ~ 3 which is three times larger than for the high-temperature uncorrelated 
fiux-fiuid but still ~ 40 times smaller than its mean field value. For g"^ = 50, (|A(G')p) has 
increased to more than half its mean field value. The insets in Fig. (|^) show the dependence 
of |A(G)p on system size for these two values of g"^. For g"^ = 30, (|A(G)p) ~ A^^^ ° as 
expected in the fiuid state while for g"^ = 50, (|A(G)p) ~ N^^'^^, consistent with the quasi- 
long-range order expected in the Abrikosov state. Fig. (0) shows that for a given system size 



(|A(G)p) increases relatively abruptly at ~ 42.5 suggesting the occurrence of a phase 
transition. 

To examine this possibility and to determine the order of the phase transition we have 
examined the dependence of the energy distribution function [|l3i p!4| on system size for 
g'^ ~ 43 and A*",^ = 80, 100, 120, 144, 168. The results are shown in Fig. @. For each systems 
size the number of Monte Carlo steps required to determine these distribution functions 
accurately exceeded 8 x 10^. For A^^ > 100 a double peak structure indicative of a first-order 
phase transition is clearly visible. For each the adjusted [l^Jl^ distribution function at 
the value of where the peaks have equal height is plotted. By extrapolating these values 
of g'^ to A*",^ = cxD as shown in the inset we estimate that a first order phase transition 
occurs at g"^ = g\j = 43.5 ± 1.0. By comparing the separations between the peak positions 
we estimate that the latent heat per flux quantum associated with the transition is ~ 
O.OlksTgl^j/ Pa ~ 0.38kBT. In Fig. (3b) we compares the (|A(G')p) distribution from values 
of the order parameter with high-energies with that from low-energies for A*",^ = 168. For 
high-energy configurations the (|A(G)p) is ~ 5.0A'^^ while for the low-energy configurations 
the distribution is peaked at ~ 0.5 demonstrating that the phase transition occurs between 
a high-energy strongly correlated vortex fluid state and a low-energy Abrikosov state. 

In Fig. (4a) we show distribution functions for /5(A) at several fixed values of Aq and 
in Fig. (4b) we show distribution functions for Aq at several fixed values of /5(A) for /3(A) 
and Aq near the values at which the phase transition takes place. The /?(A) distribution 
function is proportional to exp[{S/3{/3, Aq)/ Kb — /5A'^(7^Aq/4]. (See Eq. (pJ]). At extrema 
of the distribution dSp/djS = K^^^N^g'^Al/A.) The double peak structure apparent in the 
/5(A) distribution demonstrates that the phase transition is driven by the Sp term which 
describes the dependence of the volume in order parameter space on the degree of correlation 
in vortex positions. No similar double peak structure is seen in the Aq distribution con- 
firming that the phase transition is associated primarily with spatial correlations in vortex 
positions and hence in the superfiuid density rather than with changes in the magnitude of 
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the superconducting order parameter. 

This work was supported by the Midwest Superconductivity Consortium through D.O.E. 
grant no. DE-FG-02-90ER45427. The authors are grateful to D.P. Arovas and Kieran Mullen 
for help getting started and to S.M. Girvin for frequent useful advice. AHM acknowledges 
instructive conversations with D.A. Huse, D.R. Nelson, and Zlatko Tesanovic. 



7 



REFERENCES 
[1] A.A. Abrikosov, Zh. Eksp. Teor. Fiz. 32, 1442 (1957). 

[2] A.H. MacDonald, H. Akera, and M.R. Norman, Phys. Rev. B 45, 10147 (1992). 
[3] P. A. Lee and S.R. Shenoy, Phys. Rev. Lett. 28, 1025 (1972). 

[4] D.J. Thouless, Phys. Rev. Lett. 34, 946 (1975). G. J. Ruggeri and D. J. Thouless, J. 
Phys. F 6, 2063 (1976). 

[5] E. Brezin, A. Fujita, and S. Hikami, Phys. Rev. Lett. 65, 1949 (1990); S. Hikami, 
A. Fujita, and A. I. Larkin, Phys. Rev. B 44, 10400 (1991). 

[6] Zlatko. Tesanovic and L. Xing, Phys. Rev. Lett. 67, 2729 (1991). 

[7] Yusuke Kato and Naoto Nagaosa, Phys. Rev. B. 47, 2932 (1993). 

[8] J. A. O'Neil and M. A. Moore, Phys. Rev. Lett. 69, 2582 (1992). 

[9] Prehminary reports on these calculations have been presented previously: Jun Hu and 
A.H. MacDonald, Bull. Am. Phys. Soc. 38, 635 (1993). 

[10] If the melting transition occurs far enough from the LLL approximation becomes 
invalid and the 'London-limit' theory may apply. See for example A. Houghton, R.A. 
Pelcovits, and A. Sudbo, Phys. Rev. B 40, 6763 (1989). 

[11] For a discussion of the high-temperature perturbation expansion for XsFD{k^) see Jun 
Hu and A.H. MacDonald, submitted to Phys. Rev. B (1993). 

[12] Zlatko Tesanovic, Lei Xing, Lev Bulaevskii, Quiang Li, and M. Suenga, Phys. Rev. Lett. 
69, 3563 (1992). 

[13] J. Lee and J. M. Kosterlitz, Phys. Rev. Lett. 65, 137 (1990). 

[14] A. M. Ferrenberg and R. H. Swendsen, Phys. Rev. Lett. 61, 2635 (1988). 



8 



[15] We plan to publish further details of our results later. Jun Hu and A.H. MacDonald, in 
preparation (1993). 



9 



FIGURES 

FIG. 1. (|A(^3;)p) at = 30 (a) and = 50 (b) for a finite system with = 120. {qy = and 
T < T^^ .) The insets show the dependences of (|A(G)p) on system size at these values. The 
dashed hnes in the inset of (a) is proportional to N^^'^ while that in the inset of (b) is proportional 
to N^^'^^. {G is a member of the first shell of reciprocal lattice vectors of the Abrikosov lattice.) 
These averages were obtained from 1 ~ 2 x 10^ Monte-Carlo steps. 

FIG. 2. Dependence of {\A{q)\^) on at A^^ = 120 for g = G and for g ^ G where G is a 
reciprocal lattice vector of the Abrikosov lattice. 

FIG. 3. (a): Landau-Ginzburg energy distribution function at the finite system phase transi- 
tion point for various system sizes. Energies are in units of the mean-field condensation energy, 
N^ksTg^ / (3a- The ratio of the peak heights to the intermediate minimum grows with system 

1/2 

size but for the sizes we are able to study does not yet show the exp[cA^^' ] behavior expected at 
large N^. The inset shows the dependence of the g"^ at the phase transition on system size, (b): 
low-energy and high-energy cuts of the distribution function for |A(G)p. 

FIG. 4. Distribution functions of (a) /3[A] for several values of Aq and of (b) Aq for several 
values of /3[A]. 
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